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ABSTRACT 


The  Kidd-Stubbs  pneumatic  analogue  decompression  com¬ 
puter  is  used  for  the  real-time  prediction  of  the  safe  ascent  depth 
following  an  excursion  to  depth.  The  theoretical  model  which  is 
derived  from  the  decompression  computer  can  be  solved  using 
numerical  techniques  on  a  digital  computer,  or  can  be  solved  on 
an  analogue  computer.  This  report  describes  a  series  of  Fortran 
IV  programs  which  have  been  written  to  solve  this  theoretical 
model  under  a  variety  of  conditions.  The  basic  program  can  be 
used  to  calculate  the  safe  ascent  depth  for  a  standard  dive  or  for 
repetitive  dives,  and  if  desired,  the  safe  ascent  height  for  flying 
after  diving.  Modifications  to  the  basic  program  are  described  which 
permit  the  use  of  the  model  for  applications  such  as  the  design  of 
pneumatic  computers  and  the  analysis  of  decompression  tables 
(U) 
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INTRODUCTION 


The  Kidd-Stubhs  pneumatic  decompression  model,11'  consisting  of  four  compartments  con¬ 
nected  in  series  by  pneumatic  resistance  elements,  is  an  analogue  of  the  gas  transfer  characteristics  of  body 
tissues  during  compression  and  decompression  The  equations  of  gas  flow  between  compartments  are 
non-linear  and  cannot  be  explicitly  solved.  However,  they  are  amenable  to  solution  by  numerical  analysis 
with  a  digital  computer  or  by  recourse  to  analogue  computation  with  pneumatic  or  electronic  devices. 
The  e  techniques  have  been  described  in  an  earlier  paper,'21 


THEORY 


The  absolute  pressures  in  the  four  compartments  of  the  theoretical  model  may  be  expressed 
by  four  non-linear  first  order  dificrential  equations  describing  gas  dynamics  in  the  ‘slip-flow’  regime;'11 

dPj/dt  =  [«j  (Uj  +  I’j  +  Pj-iMPj.,  -  Pj) 

-  «j.»  <Bj.,  +fy.  +1>jH'>r  Pj,,)]  /Vj 

where  the  number  of  the  compartment  is  indicated  by  j  ( 1 ,2,3  or  4),  Pj,  Vj  are  the  absolute  pressure  and 
volume  of  the  jth  compartment,  respectively,  P()  is  the  ambient  or  driving  pressure  (input  to  the  model), 
<*j,  Uj  are  flow  constants  of  the  jth  resistor  and  Ps  =05  =  U5  =  0, 

The  constant  Oj  is  defined  by 

£*j  --  Oj  (Nd4/Lrj) 

where  N  is  the  number  of  pores  in  the  resistor  material 

L  is  the  mean  length  of  the  porous  sample  in  the  direction  of  flow 
d  is  the  mean  pore  diameter  (assuming  that  the  pores  are  of  a  circular  cross-section  ) 
and  r?  is  the  viscosity  of  the  gas 

The  value  of  Oj  may  be  determined  by  measuring  the  steady  stale  gas  flow  through  the  pneumatic 
resistor  and  by  recourse  to  the  following  equation; 

aj  =  '’is  /  taj  APs  <1Jj  +  APs  +  2Pfs) 

where  Pfs  is  the  atmospheric  pressure 

T'Aj  is  the  time  required  for  I  cm-1  of  gas  to  flow  through  the  resistor 
and  Ap  is  the  ‘driving’  pressure,  or  pressure  differential  across  the  resistor  causing  the  gas  flow. 


where 

and 


where 


and 


The  flow  constant  Bj  is  defined  by 

Bj  ~  Bj  (rj  x/T/M  /d)  (4) 

T  is  the  absolute  temperature  of  the  gas 
M  is  the  molecular  weight  of  the  gas. 

The  pressure-time  relation  for  transient  flow  of  gas  through  a  resistor-volume  combination  is 


v  npf-pi)(B  +  pf  +  rc)'l 
’  a(B  +  2Pf)  L(Pt-  Pc)(B  +  Pf  +  P))J 

a,  B  are  flow  constants  of  the  resistor, 

V  is  the  compartment  volume, 

P]  is  the  initial  compartment  pressure  at  time  t  =  0, 
Pe  is  the  compartment  pressure  at  time  t, 

P,  i:  the  final  compartment  pressure. 


(5) 


The  half-time,  T/4,  of  this  combination  is  defined  as  the  time  required  for  the  compartment  pressure  to 
decrease  from  u  value  of  P(-  to  (Pj  +  Pf)/2  and  can  be  expressed  as 

V  f2B  +  3Pf  +  Pil 

T,/l  a  (B  +  2Pf)  ln  [.  B  +  Pf  +Pj  ■"  (6) 

For  this  decompression  model,  the  half-time  is  specified  for  a  compartmental  pressure  decrease  from  50 
psig  to  25  psig.  Table  1  shows  values  of  a,  B  and  Ty,  for  several  gases  and  gas  mixtures  pertinent  to  a 
pneumatic  resistor  with  a  mean  pore  diameter  of  0.125  jam  tor  pressures  in  units  of  psi. 


TABLE  1 


Values  of  a,  B  and  T 'A  for  several  gases  and  gas  mixtures  pertinent  to  a 
pneumatic  resistor  with  a  mean  pore  diameter  of  O.I25/am. 


Gas 

M 

(amu) 

V/i 

(minutes/ 

V 

(c.c.) 

0! 

XI0~* 

B 

°2 

32 

22.4 

28.9 

4.575 

130.2 

20/80  02/N2 

28.8 

20.8 

28.9 

5.133 

122.3 

n2 

28 

20,4 

28.9 

5.309 

1 19.9 

20/80  02/lle 

9.6 

14.4 

28.9 

4.725 

230.2 

10/90  02/l  !e 

6.8 

12.5 

28.9 

4.740 

272.5 
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NUMERICAL  ANALYSIS 


Equation  (I)  consists  of  a  sc',  of  lirst-ordcr  simultaneous  differential  equations  of  the  form 
dy/dx  =  f(x,y,z) 
i//d.x  =  g<x,y,z) 

where  x  is  the  independent  variable  and  y  and  z  are  dependent  variables.  Such  a  system  is  amenable  to 
numerical  solution  by  a  Runge-Kutta  method.  The  technique  used  in  this  work  is  that  due  to  Gill^* 
characterized  by  the  following  set  of  equations,  pertinent  to  the  first  dependent  variable: 

y„.l  =yn  -E  Iky  1  +2(1  x/Wky2  +  :{l  +N/T)ky3  +  ky4]/6  +  0(h5)  ...  (7) 

where  It  -  step  size 

kyl  =  >'„> 

ky2  =hf(xn  +  ‘/4h.yn  +  !*ky|) 

ky3  =hf(xn  +  S4h.yn  +  [-  ,^  +  vr^)kyl  +[i  -  x/¥]ky2) 
ky4  =  hf  (xn  +  h,  yn  -  \/  'A  ky2  f  [  I  +  s/77 1  ky3 
This  method  is  applied  to  the  other  dependent  variables  in  a  similar  fashion. 

After  the  four  compartment  pressures  are  calculated  by  the  above  method,  the  safe  ascent 
depth  for  decompression  is  determined  from  the  greatest  of  these  pressures  divided  by  a  supersaturation 
ratio.  Experimentally,  it  has  been  found  that  the  human  body  can  withstand  some  degree  of  supersaturation 
of  inert  gas  without  experiencing  decompression  sickness.  The  supersaturation  ratio  is  defined  by  the 
equation 

kt  =  iypSA  . . .  (8) 

where  l’T  is  the  greatest  absolute  compartment  pressure 

and  l’SA  is  the  pressure  to  which  ascent  may  be  made  safely, 

Experiments  involving  more  than  4000  man-dives  at  IXTHM  to  depths  as  great  as  300  ft 
of  seawater  have  indicated  that  K  r  can  be  expressed  {for  l*|-  in  ft  of  seawater)  as 

K,- 1.385/(1  I3.7/I’r)  ...  (0) 

that  is,  as  a  function  of  the  absolute  pressure.  Table  2  shows  the  relation  between  the  supersaturation 
ratio  and  the  greatest  compartment  pressure  in  this  decompression  model. 

The  safe  ascent  depth  in  gauge  pressure,  l)SA ,  can  then  be  determined  from 


l>SA  ^  'A  •••  HO) 

where  l’A  is  the  barometric  pressure  at  the  water  surface.  The  equation  can  be  further  simplified  to 


■  (ID 
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where 

and 


R  is  a  dimensionless  decompression  factor 


Pg|,|:  is  a  pressure  offset,  Expressed  in  feet  of  seawater,  the  equation  has  the  form  (for 
PA  =  33  ft  (abs.)) 


j  -Jt. 

JsA  L385 


42.9 


. . . ( i !  a) 


TABLE  2 

PRESSURE-DEPENDENT  SUPER  SATURATION  RATiO 


Pt 

ft  seawater 

Pt-Pa 

ft  seawater 

dsa 

ft  seawater 

kt 

33 

0 

-19.07 

2.37 

59.4 

26.4 

0 

1.80 

99 

66 

28.6 

1.61 

132 

99 

52.4 

1.55 

165 

132 

76.2 

1.51 

198 

165 

100.1 

1.40 

231 

198 

123.9 

1.47 

264 

231 

147.7 

1.46 

297 

264 

171.5 

1.45 

330 

297 

195.4 

1 .445 

363 

330 

219.2 

1.44 

DIGITAL  COMPUTER  PROGRAMS 


.several  series  of  digital  computer  programs  .veto  written  in  FORTRAN  IV  for  use  on  a  PDP-9T 
computer  at  DCIFM.  All  programs  utilize  the  same  Rungc-Kiitta  calculation  method  for  determination 
of  compartment  pressures  and  the  same  safe  asccr.t  criterion. 


IDENTICAL  COMPARTMENT  PROGRAMS 


This  series  of  programs  concerns  decompression  models  in  which  all  compartments  and  pnett 
matie  resistors  are  identical.  The  input  and  output  data  files  are  similar. 


(a)  Basic  Decompression  Program  DOS: 

(I )  Options 

The  basic  program  in  this  series  is  D6S,  It  incorporates  most  of  the  features  embodied 
in  the  five  programs  described  in  Reference  2  and  operates  within  an  8K  core  of  the  PDP-9T.  Appendix 
I  contains  the  D6S  source  program  which  has  several  options  that  may  be  selected  by  a  two  digit 
integer  called  KEY  at  the  beginning  of  the  input  data  file: 

(i)  Standard  Dive:  KEY  =  00 

If  KEY  is  any  number  other  than  09  or  those  listed  below,  the  diver's  excursion  depth 
is  considered  as  a  standard  dive  in  which  descent  occurs  at  a  constant  rate,  e  g.,  60  ft/min.  to 
the  maximum  depth.  Upon  leaving  this  depth  after  a  specified  time,  the  diver  is  considered 
to  ascend  at  a  constant  rate  to  a  depth  from  which  his  ascent  is  governed  by  the  safe  ascent 
criterion  defined  in  equation  (10).  Programmed  stops  can  be  included  in  the  descent  to 
adequately  represent  a  working  diver's  excursion  to  depth. 

(ii)  Impulse  Dive:  KEY  =04 

This  type  of  dive  is  used  specifically  for  calibration  of  pneumatic  decompression 
computers  and  consists  of  a  step  impulse  of  pressure  applied  to  the  model  at  time  zero  which 
is  maintained  for  a  duration  time  I  at  which  the  pressure  is  returned  immediately  to  sea  level. 
Compartment  pressures  are  continuously  calculated  until  sea  level  is  reached. 

(iii)  Flying  after  Diving:  KEY  =  05 

This  is  a  standard  dive  in  which  calculations  of  compartment  pressures  are  continued 
after  the  model  has  ascended  to  sea  level  with  the  provisos  that  the  ambient  pressure  is 
maintained  at  sea  level  and  the  safe  ascent  'depth'  is  expressed  in  terms  of  feet  of  altitude. 
This  will  be  the  height  above  sea  level  to  which  the  model  may  safely  ascend  instantaneously 
tit  any  time  alter  surfacing. 

(iv)  Repetitive  Dive:  KEY  :  10.  20,  20,  .  .  . 

The  results  of  several  excursions  of  the  model  to  various  depths  can  he  obtained  with 
this  option  of  the  DOS  program.  It  consists  of  a  standard  dive  in  which  the  diver  ascends  to 
the  surface  via  the  safe  ascent  criterion  with  computation  of  safe  ascent  continuing  for  the 
duration  of  the  diver's  stay  at  the  surface.  Coinpuitment  pressures  at  the  etui  of  this  time  are 
used  as  initial  conditions  for  the  subsequent  dive .  The  procedure  is  repeated  N  times,  where 
N  =  KEY/ 10,  finally  finishing  at  sea  level. 
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(v)  Repetitive  Dive  Followed  by  l-lying  after  Diving:  K l-‘ Y  =  1 5,  25,  35,  .  . 

This  option  is  similar  to  that  lor  repetitive  diving  with  the  difference  that  after  the  diver 
attains  sea  level  on  the  last  excursion  to  depth,  the  safe  ascent  applicable  to  his  stay  on  the 
surface  is  expressed  in  feet  of  altitude  above  sea  level  and  is  pertinent  to  flying  after 
repetitive  diving, 

(2)  Input  Data  File  Format 

The  input  data  file  for  the  basic  decompression  program  D6S  is  denoted  as  NLDIV  SRC  and 
must  be  loaded  on  the  PDP-9T  computer  disc  for  operation  of  DbS.Tlie  following  format  must  be  used  in 
•he  composition  of  this  file: 

line  number 

I 

i 

3 

4 

5 

6 

7 

8 
9 
10 


subject 

KEY,TITL 
A,  B,  R 

DTL,  TMAX,  PTOUT,  RDES,  RASC 
T,  G 1 ,  G 
T,Gl 
-T 

TTOP 
T,  G I 
-T 

KEY,TITL 


format 

12,  1 4A5 
H8.4,  E5.I.F5.3 
5F4.I 
6F5.I 


The  nature  of  the  KEY  symbol  has  been  explained  in  the  preceding  section.  TITL  is  any 
desired  title  comprising  70  or  fewer  alphanumeric  characters.  A  is  a  flow  constant  of  the  pneumatic 
resistor  related  to  a  by  the  equation 


and  B  is  the  flow  constant  defined  earlier  in  this  paper.  R  is  the  dimensionless  computer  constant  required 
in  the  definition  of  the  supersaturation  ratio.  The  symbols  DTL  and  TMAX  denote  the  maximum  slep 
time  in  the  Runge-Kutta  computation  technique  expressed  in  minutes  and  the  maximum  length  of  the 
computer  run,  respectively.  PTOUT  is  the  printout  interval,  RDES  is  the  descent  rate  and  RASC  is  the 
ascent  rate  to  the  depth  of  safe  ascent. 

Line  4  represents  the  initial  pressure-time  stale  of  the  model  where  T  is  the  initial  time  and  Gl 
is  the  ambient  pressure  expressed  in  gauge  units  of  fee t  of  seawater.  G  represents  an  array  containing  the 
four  initial  compartment  pressures  which  need  not  be  equal  to  Gl.  If  all  are  equal  to  G  I,  then  there  is  no 
requirement  for  these  values  to  be  entered  in  the  input  file;  they  are  established  as  equal  to  Cl  I  by  the 
program  D6S. 

Line  5  represents  the  point  in  the  pressure -time  profile  at  which  the  diver  begins  Ins  ascent 
where  T  is  the  time  for  this  event  and  Gl  is  the  bottom  gauge  pressure.  The  program  D6S  calculates  the 
pressure -time  history  of  the  model  to  this  point  assuming  that  the  descent  was  made  at  a  specified 
constant  rate  and  that  the  diver  remained  at  the  bottom  gauge  pressure  specified  by  ( 1 1  in  line  5.  It 
a  stop  or  delay  is  desired  during  dcsceni  on  a  standard  dive,  the  time  that  the  diver  leaves  that  stop  and 
the  gauge  pressure  of  the  stop  can  be  inserted  immediately  before  line  5  in  the  same  format. 

When  a  negative  value  is  read  as  in  line  0,  computations  for  the  ascent  phase  of  the  dive 
commence  and  continue  until  the  diver  attains  sea  level  in  the  ease  of  the  standard  (live  or  until  (8 
hours  of  altitude  calculations  are  completed  in  the  ease  of  the  flying  after  diving  option. 
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l  ines  7  through  ‘I  are  required  only  it  the  repetitive  dive  option  is  i.-ed  T  TOP  is  the  duration 
o("  time  on  the  surface  between  the  completion  of  one  excursion  and  the  commencement  of  another.  T 
is  the  time  that  the  diver  commences  Ins  ascent  to  the  surface  from  the  hot  loin  gauge  pressure  t  I  as 
measured  from  the  start  of  the  current  excursion.  The  sequence  can  he  repeated  for  the  number  ol  dives 
required. 

Computations  pertinent  to  another  independent  dive  can  he  initiated  in  the  same  input  data 
file  by  entering  a  new  KI:Y  and  TITI  .  When  no  more  dive  profiles  are  required,  the  input  data  file  can 
be  terminated  by  setting  Kl  Y  equal  to  0*1. 

An  example  of  the  input  data  file  for  a  standard  dive  with  the  “flying  after  diving"  option  is 
shown  in  Appendix  3. 

(3)  Output  Data  f  ile  f  ormat 

The  output  of  the  basic  dccompicssion  program  is  presented  in  a  file  denoted  as  NLDIV 
LST.  The  format  of  the  presentation  consists  of  TITL.  A,  If,  R,  KASC,  RDhS  and  the  following 
columns  of  data: 

( i >  the  elapsed  time  of  the  dive. 

(ii)  the  depth  of  the  diver, 

(iii)  the  altitude  of  safe  ascent  (if  required!. 

(iv)  the  safe  ascent  depths  in  feet  of  seawater  pertinent  to  each  of  the  four  compartments. 

Appendix  3  shows  an  example  of  the  output  file. 

flic  basic  decompression  program  presents  output  on  the  teletype  or  on  the  line  printer  when 
used  in  NK  core  of  the  I’DI’-'HT  compute  system.  The  peripheral  assignments  which  must  be  made  are 

for  teletype  A  DKCI  3/T f  7. 

for  line  printer  A  DKCi  3/1.1’  7. 

(4)  Snbpiograms  Required 

One  subprogram,  fXlT,  is  required  to  operate  the  l)6S  basic  decompression  program. 


(b)  Abbreviated  Basic  Decompression  Program  DftSA: 

The  abbreviated  basic  program,  D(>SA,  permits  use  of  the  teletype,  the  tine  printer,  or  the  disk 
for  the  output  tile  when  utilizing  HK  core  of  (lie  I’DI’-OT  computer  system.  The  only  difference  from  DOS 
is  that  altitude  calculations  are  not  performed.  However,  il  the  flying  after  diving  option  is  required,  the 
safe  ascent  gauge  pressure  values  are  presented.  Consequently  the  output  data  file  differs  from  that 
lor  DCS  through  the  omission  of  an  altitude  column  but  the  input  data  files  are  identical  in  format. 
Appendix  4  shows  an  example  of  the  input  and  output  files  lor  an  impulse  dive  using  this  program. 


The  peripheral  assignments  to  he  made  to  obtain  the  required  output  mode  are  the  same  as  for 
D6S  with  the  addition  of 

disk  output  ADKB0  l/DKBI  4.  3,7. 

I  lie  abbreviated  program  DbSA  is  mot  commonly  used  because  the  living  after  diving  option  is  rarely 
required  and  because  there  is  greater  flexibility  in  output  mode. 
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(c)  Basic  Programs  with  Variable  Descent  M6S  and  M6SA: 

These  programs  are  basic  decompression  programs  similar  to  l)6S  and  D6SA  respectively  in 
which  the  rate  of  descent  can  be  varied  by  specifying  the  pressure-time  profile  at  closely  space  time  intervals 
during  the  descent  phase.  e.g. 


T 

Gl 

G.O 

0.0 

0.1 

6.0 

0.2 

12.0 

0.3 

1K.0 

0.4 

30.0 

0.5 

40.0 

0.6 

50.0 

This  type  of  descent  is  a  ‘stepwise’  approximation  to  a  continuous  pressure  change  but  if  the  intervals 
in  time  are  made  very  small,  little  error  is  made  in  the  use  of  this  approximation.  Peripheral  assignments 
for  M6S  and  M6SA  are  the  same  as  for  D6S  and  D6SA  respectively. 


(d)  Program  for  Analysis  of  Dive  Tables  D6ST. 

The  analysis  of  diving  tables  with  the  Kidd-Stubbs  decompression  model  requires  the  deter¬ 
mination  of  compartment  pressures  when  the  model  is  constrained  to  a  ‘step’  decompression  profile 
imposed  by  a  diving  table.  The  descent  and  ascent  rates  between  steps  on  the  decompression  profile 
are  established  by  RDES  and  RASC  respectively.  The  input  data  file,  NLD1V  SRC.  has  the  format: 

KEY,  TITL 
A,  B,  R 

DTL,  TMAX,  PTOUT,  RDES,  RASC 
T,  Cl,  G 

T,  Cl  (Time  at  which  pressure  maximum  is  reduced) 

T,  Cl  (Time  at  depth  of  first  stop) 

T,  G1  (Time  at  depth  of  second  stop) 


T,  Cl  (Time  at  depth  of  last  stop) 

-T 

09 

The  time  at  which  the  pressure  maximum  is  reduced  includes  the  descent  time  but  the  duration  at  each 
stop  docs  not  include  the  ascent  time  between  each  stop.  The  program  need  not  he  terminated  with 
KEY  =09  so  that  repetitive  diving  tables  may  be  analyzed.  Peripheral  assignments  are  the  same  as  for  D6SA. 

Appendix  5  shows  an  example  of  a  decompression  table  which  has  been  analyzed  with  this 
program.  If  the  safe  ascent  depth  has  not  reached  sea  level  at  the  termination  of  the  last  stop  before  surfac¬ 
ing  (an  unsafe  decompression  situation),  the  program  continues  to  calculate  the  safe  ascent  depth  until 
sea  level  is  reached  as  shown  in  the  example. 
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(e)  Modified  D6S.A  Programs: 

(I)  D6SAC2:  This  program  is  pertinent  to  operation  of  the  Kidd-Stubbs  decompression  model 
with  compensation  for  changes  in  atmospheric  pressure. The  computation  process  is  based  on  gauge  pressure 
units  and  a  How  constant  if,  defined  as 


H'  =  B  +2Pa  ...  (13) 

where  l’A  is  the  barometric  pressure. 

P,\  is  inserted  in  the  input  data  file.  NLDIV  SRC,  immediately  after  RASC.  Th1*  formula  for 
safe  ascent  depth  also  incorporates  changes  in  atmos>  heric  pressute,  viz 


I) 


SA  ~ 


-  0.578  l> 

A 


...  (14) 


where  P,  is  the  gauge  pressure  of  the  controlling  compartment. 


(2)  D6SMC.  This  program  is  a  modification  of  the  DbSA  program  that  is  applicable  to  a  Kidd- 
Stubbs  decompression  model  consisting  of  up  to  nine  identical  compartments.  The  number  of  compart¬ 
ments  desired  in  any  application  of  this  program  is  entered  into  the  input  data  file,  NLDIV  SRC.  as  a  two- 
digit  number  immediately  after  R. 

(3)  D6SAT:  A  combination  of  D6SA  and  D6ST,  this  program  is  useful  for  comparison  of  a 
standard  Kidd-Stubbs  decompression  model  dive  with  a  ‘table’  dive  for  the  same  hyperbaric  exposure 
prior  to  ascent.  It  produces  a  series  of  comparisons  for  dives  to  the  same  depth  but  for  different  time 
exposures,  e.g.  150  ft  for  10,  2C,  30  .  .  .  minutes.  The  input  data  file.  NLDIV  SRC,  has  the  following 
format  pertinent  to  standard  dives. 


line  number 


2 

3 

4 

5 

6 

7 

8 


subject  formal 

KLY.TITL 
A,  B,  R 

DTL.TMAX,  PTOUT.  RM  S.  RASC 
T.Gl.G 

T|  ,GI ,  KST  2F5.1, 12 

-T 

KEY  +  10.TITL 
T2,G1.KST 


KST  in  line  5  is  an  integer.  If  it  is  zero  then  the  program  operates  as  either  D6SA  when  KEY 
=  00  (pertaining  to  standard  dives)  or  as  D6ST  when  KEY  =  01  (pertaining  to  diving  tables).  If  KST  is  a 
positive  integer,  the  program  stores  the  computer  compartment  pressures  at  time  T( ,  for  later  use.  When 
the  dive  is  completed,  KEY  is  read  again  and  if  it  is  10  or  II.  computation  of  compartment  pressures 
resume  at  time  T( ,  using  the  stored  values,  and  continue  until  time  T2  when  the  ascent  phase  of  the  new 
jive  commences.  This  value  of  time  and  the  compartment  pressures  are  stored  by  the  computer  as 
starting  conditions  for  the  next  dive.  This  process  results  in  a  considerable  saving  of  computation  time  since 
the  compartment  pressures  pertinent  to  the  descent  phase  and  exposure  time  at  maximum  pressure  arc 
calculated  only  once. 
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Termination  of  the  series  of  dives  is  accomplished  by  setting  KEY  =0*)  and  initiation  of  a  new 
series  of  dives  occurs  when  KEY  =  00  or  01.  In  the  latter  case,  lines  2  through  4  of  the  input  data  file 
must  be  repeated.  Output  from  this  program  is  obtained  either  via  line  printer  or  via  teletype  utilizing  8K 
core  of  the  PDP-9T  computer  system. 


NON-IDENTICAL  COMPARTMENT  PROGRAMS 


In  the  construction  of  pneumatic  analogue  decompression  computers  it  is  not  possible  to 
ensure  that  all  compartments  and  pneumatic  resistors  are  identical.  A  series  of  programs  have  been  written 
to  permit  specification  of  the  exact  values  of  such  parameters  as  A,  B,  V,  and  T,^  pertinent  to  each  resistor- 
volume  configuration.  These  programs  are  intended  for  validation  of  pneumatic  analogue  decompression 
computer  performance. 

(a)  D6SVA  =  The  specification  for  each  resistor-volume  compartment  is  in  terms  of 

a,  B  and  V. 

(b)  D6SVA1  =  The  specification  of  each  resistor-volume  compartment  is  in  terms  of 

T a  ,  Ti^  and  B. 

(c)  D6SVA2  =  The  specification  of  each  resistor-volume  compartment  is  in  terms  of 

V,  T^  and  B. 

The  input  data  file  format  for  D6SVA  is  as  follows: 


line  number  subject  format 


1 

KEY,TITL 

12,  14As 

2 

A(  1 ),  A(2),  A(3),  A(4) 

4E8.4 

3 

B(l),  B(2),  B(3),  B(4) 

4F5.I 

4 

V(l),  V(2),  V(3),  V(4) 

4F5.2 

5 

R 

F5.3 

6 

DTL,TMAX,  PTOUT,  RDES,  RASC 

5F5.I 

7 

TGI  ,C 

6F5.1 

8 

T,GI 

9 

T 

10 

KEY,  TITL 

F  or  the  other  two  programs  in  this  series,  the  input  data  format  is  the  same  with  the  appropriate 
changes  in  lines  2  through  4. 


ON-LINE  PROGRAM  AT  HYPERBARIC  CHAMBER:  DIVER 


I  I 


The  basic  decompression  program  concerning  four  identical  compartments  has  been  adapted 
for  operation  in  real-time  from  the  hyperbaric  chamber  and  is  shown  in  Appendix  2.  The  output 
voltage  from  a  pressure  transducer  in  the  300  ft  hyperbaric  chamber  at  DCIEM  is  transmitted  as  real-time 
digital  information  into  the  l’l)P-‘)T  computer  system  via  an  analogue-to-digital  converter.  In  operation 
of  the  DIVER  program,  pressure  in  the  chamber  is  sampled  every  10  seconds  and  calculation  of  compart¬ 
ment  pressures  is  made  with  the  Runge-Kutta  routine. 

A  conversational  mode  is  used  to  start  this  pro  -  rai  >,  to  input  data  via  the  teletype  and  to 
instruct  the  chantbci  computer  operator  as  to  the  --peraiiv  >  of  the  pr  'g.ao  .  Output  of  the  elapsed  time 
in  the  dive,  pressure  in  the  chamber,  safe  ascent  depth,  the  number  of  the  controlling  compartment  and 
the  pressures  in  each  of  the  four  compartments  occurs  every  minute  via  teletype.  If  the  pressure  in  the 
hyperbaric  chamber  is  permitted  to  drop  below  the  predicted  safe  ascent  depth,  the  program  prints  out  the 
word  ‘VIOLATION’  on  th;  teletype  as  a  warning  to  the  chamber  controller.  At  the  end  of  a  dive  the 
sampling  of  the  pressure  may  be  terminated  by  typing  CONTROL  P.  Theoretical  calculations  are  continued 
in  ■•'-.mputer  time  representative  of  12  hours  with  printout  corresponding  to  one  hour  intervals. 

Two  related  versions  of  DIVER  have  been  written.  DIVER1  presents  only  such  basic  infor¬ 
mation  n  output  to  the  chamber  controller  as  clasped  time,  chamber  pressure,  safe  ascent  depth  and 
the  number  of  the  controlling  compartment.  DIVER2  is  a  version  that  presents  the  output  of  D!  /ERI 
on  teletype  and  stores  the  rcal-timc  output  from  L-iVER  on  DECTAI’E. 

In  addition  to  the  subprogram,  EXIT,  these  programs  require  a  monitor  program,  ADF4TS, 
which  allows  the  use  of  the  analogue-to-digital  converter. 
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APPENDIX  1 

BASIC  DECOMPRESSION  PROGRAM  D6S 
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C  FILE  D6S  *  BASIC  6S  MODEL  WITH  RAMP  DESCENT 

C  (WITH  ALTITUDE  CALCULATIONS) 

C 

REAL  INDATA  (2 ) , OUTDAT (2 ) 

LOGICAL  FOUND, C.KTEST 

COMMON  TITL(14),P(4),G(4),GBIG(4),AK(4,4),F(5) 

DATA  INDATAC 1 ) , I NDATA ( 2) /5HNLDI V, 4H  SRC/ 

DATA  OUTDATC ’.  )  ,OUTDAT(2)/5HNLDIV,4H  LST/ 

DATA  Z 1,Z2,Z3/. 207136781 , .292893219, .7071 06781 / 

DATA  Z4,Z5,Z6/1 .707106781 ,  .585736433,3.414213562/ 

C 

C  OPEN  DATA  FILE 
C 

CALL  SEEK  (3, INDATA) 

CALL  DLETE  ( 7 , OUTDAT , FOU ND) 

CALL  ENTER  (7, OUTDAT) 

1  READ  (3, 101 )  KEY,  TITL 

101  FORMAT  (12,  14A5) 

IF  (KEY. ME. 9)  GO  TO  2 
CALL  CLOSE  (3) 

CALL  CLOSE  (7) 

CALL  EXIT 
C 

2  WRITE  (7,102)  TITL 

102  FORMAT  ( 1H1 , 14A5) 

READ( 3 , 103)  A , 9  , R 

103  FORMAT  ( E3.4,F5.1,F5.3) 

READ(3,1 12)  DTL,TMAX,PTOUT,RDES,RASC 
112  F0RMAT(5F4.1 ) 

WRITE  (7, 103)  A,  B,  R,  DTL 

108  FORMAT(  1  H3 , 4X4HA  =  ,F9.B,5X4HB  r  ,F5.1,5X4HR  =  ,F5.3/ 

1  1 H0 , 3X1 6HMAX  STEP  TIME  =  ,Fa.2,4H  MIN) 

IF  (KEY.EQ.4)  GO  TO  125 
WRITE  (7,  120)  RASC ,RDES 

120  FORMAT (4 X5 HA5C  =F5.1,7H  FT/MI N6X6HDESC  =F5.1,7H  FT/MI N/) 

125  KTEST=5*(KEY/5) .EQ.KEY 

TMAXUTMAX 
WRITE  (7,109) 

109  FORMAT( lH0,3i  H  TOTAL  DIVE  ACTUAL  ALTITUDE 

1  6X35  HSAFE  ASCENT  DEPTHS  (FT  OF  SEA  WATER)  /2X1 0HTIME  (MIN) 

2  2 1 H  DEPTH  (FT)  (1000  FT)9X1HI , 8X 1 H2 , 8X 1 H3 , 8X1 H4/) 
DTUDTL*.  1 

C 

C  INITIAL  TIME  AND  PRESSURES 
C 

6  READ  (3,134)  T,  61, G 

104  FORMAT  '6F5.1) 

K  =  1 

FOUND=G( I ).EQ.0. 

DO  7  1=1,4 
IF  (FOUND)  G  ( I )  r  G 1 
P(I)  =  G ( I )  +  33. 

WrGBIG(I)  =  P(I)/R-42.9 
IF  (GBIG(K).LT.W)  K=I 


7 


14 


PI  -  G1  +  33. 
PBIGK=GBIG(K)  +  33. 

TT  :  T 
C=. FALSE. 

NS  3  =  5 
GO  TO  30 
C 

25  PMIN  =  0. 

KPT  =  0 
C 

C  READ  IN  STEP  TIME  AND  PRESSURES 
C 

11  DT  =  DTL 

PLAST  =  PI 
READ  (3,104)  T,  G1 
17  (T.LT.0.)  GO  TO  12 
IF (NS2.EQ.2)  T  =  T+TT 
NS2~  1 

TL-T  -  .01 

IF  (KEY.EQ.4)  GO  TO  59 
DP  =  . l+RDES 


W2-';l  +  33.  - 

DP 

56 

IF  (Pi.LT.W2) 

GO  TO  58 

59 

PUGH-33. 

NS  =  5 

GO  TO  57 

58 

DT=0.1 

NS:  1 

PUPl+DP 

C 

C  TEST  FOR  DESCENT  PRESSURE  AND  TIME  CHANGES 
C 

57  IF  (PI. NE. PLAST)  DT=DT1 

19  IF  (TT.LT.TL)  GO  TO  10 

IF((TTEST.LT..05) .UR.(TTEST.GT.(PTOUT-.05)))GO  TO  11 
43  NS3=1 

GO  TO  30 
C 

C  PART  OF  RUNGE-KUTTA  ROUTINE 
C 


22 

DO  14  I  :  1,4 

W:GSIG(I)  :  P(I)/R- 

14 

IF  (GBIG(K)  .LT.W)  K 

PBIGK=GBIG(K)+  33. 

TT  =  TT  +  DT 

TTEST  =  AMOD  (TT,PTOUT) 

.5  *DT 

IF((TTEST.GT.W). AND. (TTEST. LE .( PTOUT-W) ) )  GO  TO  32 
IF  (KPT.EQ.l)  DT=DTL 

MC  1-0 

77  GO  TO  30 

32  IFUPBIGK.  GT.  PMIN)  .AND.  (TT.LT.TMAX1  ))GO  TO  ( 56 , 1 0 , 1  8 , 1  7, 57) ,  NS 

NS3=3 
GO  TO  30 


c 

C  TESTS  FOR  REPETITIVE  AND  FLYING  AFTER  DIVING 
C 

33  IF  (KEY.GE.10)  GO  TO  (23,26),NS2 

IF  (KEY.EQ.5)  GO  TO  20 
GO  TO  I 
C 

12  PMIN  =  33. 

C 

C  TEST  FOR  IMPULSE  ASCENT 
C 

IF  (KEY.NE.4)  GO  TO  13 
C 

C  IMPULSE  ASCENT 
C 

21  NS=2 

24  Pl=33. 

C  =  .FALSE. 

GO  TO  10 
C 

C  STANDARD  ASCENT  TO  FIRST  STOP 
C 


13 

DT 

=  0. 

1 

DP 

=  0. 

1  *RASC 

NS: 

:3 

18 

PI 

=  PI 

-DP 

IF 

(PI  . 

GT. PBIGK)  GO  TO  10 

C 

C  CONTINUOUS  CONTROLLED  ASCENT 
C 

C  =  .TRUE. 

KPT  =  1 
NS=NS3  =  4 
PI  =  P8IGK 
GO  TO  30 

17  PI  =  PBIGK 

GO  TO  10 
C 

C  FLYING  AFTFR  DIVING  ROUTINE 
C 

20  KEY  =  0 

TMAXUTT+1080. 

PTOU  T= 1 5  .0 
GO  TO  41 
C 

C  REPETITIVE  DIVING  ROUTINE  -  SURFACE  INTERVAL 


C 

W  IS  TIME  AT  SURFACE 

C 

23 

READ  (3,105)  W 

105 

FORMAT  (F5.I) 

TMAX1  =  TT  +  W 

NS2-2 

4! 

PMI N=0 . 

GO  TO  21 

C 

C  DESCENT  AND  DIVE  REPETITION 
C 

26  TMAX1  =  TMAX 

KEY  =  KEY  -  10 

GO  TO  25 
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C 

C  OUTPUT  PRINT  ROUTINE! 

C 

30  WrPBIGK 

IF  ((KEY.EQ.4).OR.(T.GE.0.).OR.(  .NOT  .C) )  W=P1 
GB  =  W-33. 

38  WRITE  (7,110)  TT,GB,K,(GBIG(I),I=l,4) 

1 10  FORMAT (3XF7,2,5XF6.2, 1 4X1 1,4(3XF6,2)) 

IF  (.NOT.KTEST)  GO  TO  44 
IF  ((GBIG(K)+.0001 ).GE.0.)  GO  TO  44 
W=1 45.5300* <1 ,-(PBIGK/33.)**.l903) 

WRITE  (7,115)  W 

115  FORMAT  ( 1H+,25X,F5.2) 

44  GO  TO  (11, 32, 33, 10, 25), NS3 

C 

C  RUNGE-KUTTA  ROUTINE 
C 

10  J  =  1 

F(5)=0. 

X=P1 

DO  401  1=1,4 
Y=P(I) 

GO  TO  450 

401  CONTINUE 
GO  TO  500 

402  IF  (C)  X=X+AK(l ,K)*,5 
DO  403  1=1,4 
Y=P(I)+AK(1 ,I)*.5 

GO  TO  450 

403  CONTINUE 
GO  TO  500 

404  IF  (C)  X=X+Z1*AK(1 ,K)+Z2*AK(2,K) 

DO  405  1=1 ,4 

Y=P(I)+Zl*AK(l,I)+Z2*AK(2,I) 

GO  TO  450 

405  CONTINUE 
GO  TO  500 

406  IF  (C)  X=X-Z3*AK(2,K)+Z4*AK(3,K) 

DO  437  I = l ,4 

Y=P( I)-Z3*AK (2, 1 )+Z4*AK (3 , 1 ) 

GO  TO  450 

407  CONTINUE 
GO  TO  500 

408  DO  409  1=1,4 

409  P(I)=P(I)+(AK(l, I)+Z5*AK(2,I)+Z6*AK(3,I)+ 

1  AK(4,I))/6. 

GO  TO  22 

450  F(I)=A*(X-Y)*(B+X+Y) 

X=Y 

GO  TO  (401 ,403,405,407)  ,J 

500  DO  501  1=1,4 

501  AK(J,I)=(F(I)-F(I+l))*DT 
X  =  P1 

J  =  J+1 

GO  TO  ( 10, 402, 404, 406, 408), J 
END 


APPENDIX  2 
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FILE  DIVER,  FOR  ONLINE,  REALTIME  OPERATION  AT  HYPERBARIC 
CHAMBER( OPERATOR  VERSION  WITH  TELETYPE  OUTPUT  ONLY) 

ALL  DATA  OUTPUT  ON  TELETYPE 

INTEGER  AHG1 ,ARG2,ARG3 

COMMON  P(4) , GBI G( 4 ) ,AK (4 , 4) , F(5 ) , TI TL(8 ) 

DATA  Z1,Z2,Z3/. 207106 781 ,  .2928932 1 9 , ,7071 06781 / 

DATA  Z 4, Z5,Z6/1 .707106781 , .585786438,3 .414213562/ 

DATA  A1 ,B1,R/. 0000791 2, 274. 5,1  .385/ 

DATA  A2, B2/. 00007281 »516. 7/ 

INSTRUCTIONS  FOR  OPERATOR 

WRITE( 7,1) 

FORMAT (/I X43HTYPE  DIVE  NO.  AS  DXXXX,  DATE,  BREATHING  GAS) 
READ(4,8) TITL 
FORMAT (8A5) 

WRI TE ( 7 , 4 ) 

FORMAT(/35H  GAS  MIX.  TYPE  1  IF  AIR,  2  IF  02-HE) 

READ(4,5)KEY 

FORMAT (II) 

WR I  TE(  7, 6) 

F0RMAT(/34H  TYPE  CHANNEL  NO.  AS  1,  2,  3,  OR  4) 

READ(4,5) ARG1 

SET  INITIAL  CONDITIONS 

IF(KEY.EQ.2)  GO  TO  90 
A=  A 1 
B  =  B1 

GO  TO  95 
A:A2 
B  =  B2 

WRITE(7, 18) 

F0RMAT(/13H  BRING  A/D  ON  LINE) 

ARG3= 1 000 

CALL  ADINIT(ARGl) 

CALL  ADREADCARGl ,ARG2,ARG3) 

CALL  ADSTOP(ARGl) 

PA=FLOAT(AR G2)/2. 

WRITEC7,  19) 

FORMAT(/44H  TURN  A/D  OFF-LINE,  PRESS  RETURN  ON  TELETYPE/) 
READ<4,28) 

FORMATO 

WRITE<7,9)TITL 

F0RMAT(//1X8A5/) 

WRITE(7,  13 )  PA 

FORMAT( 5X13HATM .  PRESS.  =F5.1,9H  FT (ABS . ) ) 

WRITE  (7,108)  A,  B,  R 

FORMAT( 5X4HA  =  F9.8,5X4HB  =  F5.1,5X4HR  =  F5.3) 
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WRITE  (7,109) 

109  F0RMAK//31H  TOTAL  DIVE  ACTUAL  SAFE  ASC 
1  6X35HC0MPARTMENT  DEPTHS  (FT  OF  SEAWATER)) 

WRITE(7, 112) 

112  FORMAT (32H  TIME  (MIN)  DEPTH(FT)  DEPTH(FT) 

1  9X1H1 , 8X1 H2 , 8X1 H3 , 8X1 H4/) 

C 

97  PTOUT=1.0 

PI  =  PA 
K  =  1 

DO  7  1=1,4 
P( I ) =PA 

7  GBIG(I)=-19.07 

NS  1=  NS2=  NS3=  1 
TT=0. 

60  TO  30 
C 

11  WRITE(7,29> 

29  FORMAT (IX, 10HSTART  DIVE) 

27  CALL  ADINIT(ARGl) 

C 

C  START  DIVE 

C 

CALL  CNTRLP 
GO  TO  33 
C 

12  CALL  ADREAD(ARG1,ARG2,ARG3) 

PI = FLOAT ( ARG2) /2 . 

dt=  float ( arg3 ) /6000 , 

GO  TO  10 
C 

C  SELECTION  OF  DEEPEST  COMPARTMENT  AND  TEST  FOR  PRINT 

C 

22  DO  14  I  =  1,4 

W=GBIG(I)=(P(I)-PA-2S.4)/R 

14  IF  (GBIG(K).LT.W)  K  =  I 
TT  =  TT  +  DT 

TTEST  =  AMOD  (TT,PTOUT) 

W=.5*DT 

IF((TTEST.GT.W). AND. (TTEST. LE . ( PTOUT-W) ) )  GO  TO  ( 1 2, 10 ,58) , NS2 
GO  TO  (1 5, 48), NS l 

15  NS3=2 

46  GO  TO  30 

C 

48  T=TT/60. 

GO  TO  49 
C 

C  OUTPUT  PRINT  ROUTINE 

C 

ir‘  T=TT 

m9  GK=GBIG(K) 

W=P1 -PA 

42  WRITE  (7,110)  T,W,GK,K,(GBIG(I), 1=1,4) 

110  F0RMAT(2X,F6.1,5XF6.1,4XF6. 1,4X12, 4(3XF6.D) 

IF(W.LT.GK)  GO  TO  31 
GO  TO  (11,12, 10,37,25), NS3 
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C 

31  WRITEC7, 16) 

36  FORMAT ( 1 X9HVI0LATI0N/) 

GO  TO  44 
C 

C  RUNGE-KUTTA  ROUTINE 

C 

10  J=1 

F(5)=0. 

X=P1 

DO  401  1=1,4 
Y=P(I) 

GO  TO  450 

431  CONTINUE 

GO  TO  500 

402  DO  433  1=1,4 
Y=P(I)+AK(1 ,I>*.5 
GO  TO  450 

403  CONTINUE 
GO  TO  503 

434  DO  405  1=1,4 

Y=P(  I)+Zt*AK(!,I)+Z2*AK(2,I) 

GO  TO  450 

405  CONTINUE 
GO  TO  500 

406  DO  407  1=1,4 
Y=P(I)-Z3*AK(2,I)+Z4*AK(3,I) 

GO  TO  453 

407  CONTINUE 
GO  TO  500 

408  DO  409  1=1 ,4 

409  P( I)=P( I)+(AK ( 1 , I)+Z5*AK(2,I>+Z6*AK(3,I)+AK(4, I) )/6. 
GO  TO  22 

450  FCI)=A*(X-Y)*(B+X+Y) 

X=Y 

GO  TO  (401 ,403, 405, 407), J 

500  DO  501  1=1,4 

501  AK(J, I)=(F( I)-F( 1+1 ))*DT 
X=P1 

J = J+ 1 

GO  TO  (10, 402, 404, 406, 408), J 
C 

c  END  OF  DIVE 

C 

33  CALL  ADSTOP  (ARG1) 

NS1=NS2=2 
NS3  =  5 
PA  =  P1 

U  TTE( 7,17) 

CA  .  CNTRLP 
GO  - o  35 
GO  TO  30 
C 
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C  THEORETICAL  CALCULATIONS  FOR  12  HOURS  AFTER  DIVE 

C 

25  WR  ITE(  7, 1$ ) PA 

WRITE(7,40) 

40  FORMAT ( 1 2H  END  OF  DIVE) 

NS3  =  4 

WRITE(7f 41) 

41  FORMAT ( 9 X5 HHOURS ) 

GO  TO  37 

C 

35  WR  ITE  ( 7 ,  17) 

17  FORMAT (/) 

GO  TO  97 
C 

37  DT= 1 .0 
PTOUT=60. 

NS2=NS3=3 

TMAX=TT+720. 

GO  TO  10 

C 

38  IF(TT.GT.TMAX)  CALL  EXIT 
GO  TO  10 

END 
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APPENDIX  3 

EXAMPLE  OF  INPUT  AND  OUTPUT  FILES  FOR  D6S  PROGRAM 
INPUT  DATA  FILE 

05  STANDARD  DIVE  200  FT  FOR  |3  MINUTES  (FLYING  OPTION) 

. 791 2E- 40 274501 33 5 
30 1 39 99900280 60 00 600 
3030000000 
3312332300 
-1 
0? 

OUTPUT  FILE 

STANDARD  DIVE  230  FT  FOR  10  MINUTES  (FLYING  OPTION) 

A  :  .00007912  B  :  274.5  R  =  1.33^ 


MAX  STEP 

TIME  =  1 .00 

MIN 

ASC  -  60 

.0  FT/MIN 

DESC  = 

60. 

0  FT/MIN 

TOTAL  DIVE 

ACTUAL 

ALTITUDE 

SAFE  ASCENT  DEPTHS 

(FT  OF 

S' EA  WATER) 

TIME  (MIN) 

DEPTH  (FT) 

(1000  FT) 

1 

2 

3 

4 

0.30 

0.00 

22.03 

1 

-19.07 

-19.07 

-19.07 

-19.07 

2.00 

120.00 

17.45 

1 

-16.14 

-19.02 

-19.07 

-19.07 

4  .00 

200.00 

5.93 

1 

-6.53 

-13.64 

-19.06 

-  19.07 

6.03 

200 .00 

1 

3.30 

-17.70 

-19  .02 

-19.07 

3.00 

200.00 

1 

12.96 

-15.26 

-18.92 

-19.07 

10.03 

200.00 

1 

21,04 

-14.42 

-18.74 

-19.05 

12.00 

80.00 

1 

23.65 

-12.36 

-18.47 

-19.03 

12.90 

22.44 

1 

22.44 

-11 .50 

-18.33 

-19.02 

14.00 

20.33 

1 

20.33 

-10.57 

-18.13 

-13.99 

16.00 

16.96 

1 

16.96 

-9.20 

-17.73 

-18.94 

18.03 

13.92 

1 

13.92 

-8.17 

-17.31 

-13.36 

20.00 

11  .21 

1 

1  1  .21 

-7.41 

-1 5.87 

-13.77 

22.00 

8.77 

1 

8.77 

-6.88 

-16.44 

-13.66 

24.00 

6.53 

1 

6.53 

-6.53 

-16.02 

-18.53 

26.00 

4.58 

1 

4.53 

-6.32 

-15,63 

-13.33 

28.03 

2.77 

1 

2.77 

-6.24 

-15.86 

-18.23 

30.00 

1.11 

1 

1  .1  1 

-6.25 

-14.92 

-13.06 

32.00 

-0.41 

0.35 

1 

-0.41 

-6.33 

-14.60 

-17.89 

32.00 

-0.41 

0.35 

1 

-0.41 

-6.33 

-14.60 

-17.3? 

45.03 

0 .00 

6.54 

1 

-7.08 

-7.79 

-13.27 

-16.65 

60  .00 

0 .30 

9.41 

2 

-11.01 

-9.77 

-12.75 

-15.37 

75.00 

0.03 

1  1  .24 

2 

-13.19 

-1 1 .37 

-12.75 

-14.46 

90.00 

0.03 

12.66 

2 

-14.51 

-12.54 

-12.95 

-13. si 

105.00 

0.30 

13.43 

i 

-15.37 

-13.40 

-  13.20 

-13. 62 

120.00 

0 . 00 

13.32 

3 

-15.95 

-14.04 

-13.47 

-13.53 

135.00 

0.00 

13.93 

4 

-16.36 

-1 4.53 

-  13.73 

-13.55 

150  .03 

0 .00 

14.07 

4 

-16.66 

-14.91 

-13.99 

-  13.66 

165.00 

0.00 

14.28 

4 

- i 6 .39 

-15.22 

-14.23 

-13.  ;2 

180.03 

0.00 

14.52 

4 

-17.07 

-15.47 

-  1  4,46 

-14,00 

195.33 

0.00 

14.78 

4 

-17.21 

-15.69 

-14.69 

-14. 20 

213.00 

0.00 

15.05 

4 

-17.33 

-15  .39 

-14.30 

-14.41 

225.30 

0.00 

15.33 

4 

-17.44 

-16.06 

-15.10 

-14.61 

\ 


APPENDIX  4 

EXAMPLE  OF  INPUT  AND  OUTPUT  FILES  FOR  D6SA  PROGRAM 
I  ".'PUT  DATA  FILE 

04  IMPULSE  DIVE  200  FT  FOR  30  MINUTES 
.73 1 2£- 4fl 27 450 1 335 
CO  1 09999335005000600 
0003003330 
0030002330 
-1 
09 

OUTPUT  FILE 

IMPULSE  DIVE  200  FT  FOR  30  MINUTES 

A  =  .00037912  3  =  274.5  R  =  1.385 


MAX  STEP  TIME  =  1 .00  MIN 
TOTAL  DIVE  ACTUAL 

CRITICAL 

SAFE  ASCENT  DEPTHS 

(FT  OF 

SEAWATER) 

TIME  (MIN) 

DEPTH  (FT) 

VOLUME 

1 

2 

3 

4 

3  ,00 

0  .00 

1 

-19.07 

-19.07 

-19.07 

-19  .07 

5.00 

200.00 

1 

7.31 

-17.20 

-18.99 

-19.37 

10.00 

200.00 

1 

27.63 

-12.47 

-13.51 

-19  .04 

15.00 

200.00 

1 

42.00 

-6.1  1 

-17.42 

-13.91 

20.00 

200.00 

1 

52.44 

0.39 

-15. 69 

-13.61 

25.00 

200.00 

1 

60.17 

7.88 

-13.36 

-13.09 

30.00 

200.00 

1 

66.03 

14.49 

-10.55 

-17.29 

35.00 

0.00 

1 

45.53 

17.89 

-7.55 

-16.20 

40.00 

3.00 

1 

32.16 

17.73 

-4.96 

-14.85 

45.00 

J  .00 

! 

22.90 

16.20 

-3.03 

-13.37 

50.00 

0.00 

1 

16.17 

14.09 

-1.70 

-11  .85 

55.00 

0.00 

2 

11.10 

1  1  .88 

-0.87 

-10.39 

60.00 

0.00 

2 

7.17 

9.75 

-0.39 

-9.02 

65.30 

0.00 

2 

4.04 

7.78 

-0.17 

-7.80 

70.00 

0.00 

2 

1  .51 

6.00 

-0.13 

-6.72 

75.00 

0.00 

2 

-0.57 

4.40 

-0.21 

-5.79 

80.03 

0.00 

2 

-2.31 

2.97 

-0.37 

-5.02 

85.00 

0.00 

2 

-3.77 

1  ,70 

-0.57 

-4.37 

90.00 

0.00 

2 

-5.01 

0.56 

-0.8  1 

-3.35 

93  .03 

0 .00 

2 

-5.67 

-0.05 

-0.97 

-3.59 

APPENDIX  5 

EXAMPLE  OF  INPUT  AND  OUTPUT  FILES  FOR  D6ST  PROGRAM 


INPUT  DATA  FILE 

0!  US  ft  DIVE  TABLE  150  FT  FOR  30  MINUTES 
. 791 2E -40271501385 
03 1 099 99 0020 060 03 600 
000303 3003 
3330001503 
33C3033230 
0324300100 
-1 
09 

OUTPUT  FILE 

US N  DIVE  TABLE  150  FT  FOR  30  MINUTES 

A  :  .00007912  B  =  274.50  R  =  1.385 
MAX  STEP  TIME  =  1.00  MIN 
4SC  =  60.0  FT/MI N  DESC  r  60.0  FT/MI N 


TOTAL  DIVE 

ACTUAL 

CRITICAL 

SAFE  ASCENT  DEPTHS 

(FT  OF 

SEAWATER) 

TIME  (MIN) 

DEPTH  (FT) 

VOLUME 

1 

2 

3 

4 

0.00 

0.00 

1 

-19.07 

-19.07 

-19.07 

-19.07 

2.00 

120.00 

1 

-16.14 

-19.02 

-19.07 

-19.07 

4.00 

150.00 

1 

-3.67 

-18.67 

-19.06 

-19.07 

6.00 

150.00 

1 

-1  .31 

-17.96 

-19.03 

-19.07 

8.00 

150.00 

1 

4.30 

-16.94 

-18.95 

-19.07 

10.00 

150.00 

1 

9.74 

-15.67 

-13.82 

-19.06 

12.03 

150.00 

1 

14.57 

-14.20 

-13.63 

-19.04 

14.00 

150.00 

1 

13.87 

-12.53 

-18.37 

-19.01 

16.00 

150.00 

1 

22.63 

-10.86 

-18.05 

-13.97 

13.00 

150.00 

1 

26.03 

-9.05 

-17.67 

-13.91 

20.00 

150.00 

1 

29.12 

-7.20 

-17.22 

-18.83 

22.00 

150.00 

1 

31  .83 

-5.33 

-16.71 

-18.74 

24.00 

150.00 

1 

34.26 

-3.46 

-16.13 

-13.61 

26.00 

150.00 

1 

36.45 

-1  .60 

-15.51 

-13.46 

28.00 

150.00 

1 

33.43 

0.23 

-14.33 

-13.29 

30.00 

150.03 

1 

40.23 

2.02 

-14.11 

-13.03 

32.00 

30.00 

l 

33.05 

3.66 

-13.35 

-17.35 

34.00 

20.00 

1 

33.03 

4  .84 

-12.57 

-  1  7.59 

36.00 

20.00 

1 

23.36 

5.53 

-11.31 

-17.29 

33.00 

20.00 

1 

25.26 

6.00 

-11.07 

-16.97 

40.00 

20.00 

22  J  7 

6.13 

-10.39 

-16.63 

40.20 

20.00 

1 

21  .89 

6.19 

-10.32 

-  1 6.60 

42.00 

10. P'1 

1 

19.15 

6.17 

-9.75 

-16.27 

44.00 

10.00 

1 

16.43 

6.00 

-9.17 

-15,90 

46.00 

13.00 

1 

14.14 

5.72 

-8.65 

-15.51 

43.00 

10.00 

1 

12.03 

5.36 

-8.13 

-15.12 

50.00 

10.00 

1 

10.26 

4.95 

-7.77 

-14.72 

52.00 

10.00 

1 

3.64 

4.50 

-7.41 

-14.32 

54.00 

10.00 

1 

7.19 

4.02 

-7.09 

-13.92 

56.00 

10  .00 

1 

5.39 

3.54 

-6.82 

-13.53 

53.00 

10.00 

1 

4.72 

3.05 

-6.59 

-13.15 

60  .00 

10.00 

1 

3.66 

2.57 

-6.40 

-  12.77 

62.00 

13.00 

1 

2.70 

2.09 

-6.23 

-12.41 

64 . 00 

10.00 

1 

1  .83 

1  .62 

-6.09 

-12.06 

64.40 

10.00 

1 

1  ,66 

1  .53 

-6.07 

-1 1 .99 

66 . 40 

0  .00 

2 

0.50 

1  .07 

-5.96 

-11.65 

63.40 

0.00 

2 

-0.55 

0.60 

-5.87 

-11.32 

70 . 43 

0.00 

2 

-1.51 

0.14 

-5.81 

-11.01 

71.40 

0.00 

2 

-1  .9. 
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-5.78 

-10.36 

